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ABSTRACT 

Modified Newtonian Dynamics (MOND) has been shown to be able to fit spiral 
galaxy rotation curves as well as giving a theoretical foundation for empirically de- 
termined scaling relations, such as the Tully- Fisher law, without the need for a dark 
matter halo. As a complementary analysis, one should investigate whether MOND 
can also reproduce the dynamics of early -type galaxies (ETGs) without dark matter. 
As a first step, we here show that MOND can indeed fit the observed central veloc- 
ity dispersion do of a large sample of ETGs assuming a simple MOND interpolating 
functions and constant anisotropy. We also show that, under some assumptions on the 
luminosity dependence of the Sersic n parameter and the stellar M/L ratio, MOND 
predicts a fundamental plane for ETGs : a log - linear relation among the effective ra- 
dius R e ff, Co and the mean effective intensity (I e )- However, we predict a tilt between 
the observed and the MOND fundamental planes. 
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1 INTRODUCTION 

The regularity of their photometric properties and the ex- 
istence of remarkable scaling relations among their observ- 
able quantities may naively suggest that early - type galaxies 
(ETGs) are well understood systems. On the contrary, un- 
derstanding their mass content and density profile is still 
a hotly debated issues mainly because of both theoretical 
shortcomings and observational difficulties. Indeed, on the 
one hand, the lack of a reliable mass tracer makes it difficult 
to constrain the gravitational potential in the outer (sup- 
posedly dark matter dominated) regions, although planetary 
nebulae surveys (Napolitano et al. 2001, 2002; Romanowsky 
et al. 2003) are trying to address this problem. On the other 
hand, the availability of a higher number of possible tracers 
in the inner regions has not improved so much the situation 
with the same data be reproducible in terms of different 
(and somewhat contrasting) scenarios because of the uncer- 
tainties on the stellar initial mass function (IMF). 

In the classical Newtonian framework, dark matter 
(DM) dominates the outer galaxy regions (see, e.g., van den 
Bosch et al. 2003a, 2003b, 2007), but weighting its contri- 
bution in the inner regions where stellar matter plays a not 



negligible role is quite difficult. As pointed out by Mamon & 
Lokas (2005a, b) , the observational data claim for a dominant 
stellar component at a radius <C -Re//, but such a conclusion 
heavily relies on the IMF choice, with the Salpeter (1995) 
and Chabrier (2001) as leading but not unique candidates. 
On large scales, the DM content has been found to be a 
strong function of both luminosity and mass (Benson et al. 
2000; Marinoni & Hudson 2002; van den Bosch et al. 2007) 
with a different behaviour between faint and bright systems. 
Looking for a similar result for the DM content at ~ Reff 
is quite controversial. On one hand, some authors (Gerhard 
et al. 2001; Borriello et al. 2003) argue for no dependence. 
In contrast, other works (Napolitano et al. 2005; Cappel- 
lari et al. 2006; Tortora et al. 2009) do find that brighter 
galaxies have a larger DM content, while a flattening and 
a possible inversion of this trend for lower mass systems, 
similar to the one observed for late -type galaxies (Persic 
et al. 1993), is still under analysis (Napolitano et al. 2005; 
Tortora et al. 2009) with no conclusive result yet obtained. 
It is worth stressing that all these results rely on dynami- 
cal analysis, i.e. they are obtained by fitting a given model 
to the observed velocity dispersion data. As such, they are 
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plagued by the mass - anisotropy degeneracy introducing an 
unknown systematic bias which is difficult to quantify. As 
an alternative, one can rely on gravitational lensing which 
directly probes the full mass content projected along the 
line of sight. In the strong lensing regime, the formation of 
Einstein rings allows to constrain the total mass projected 
within the Einstein ring independently of the model. If cou- 
pled to a measurement of the central velocity dispersion, this 
method gives interesting constraints on the mass profile in 
the inner regions. On one hand, dark matter mass fractions 
and their scaling with stellar mass and luminosity can be 
constrained assuming a model for the halo (Koopmans & 
Treu 2003; Grillo et al. 2009; Auger et al. 2009; Auger et al. 
2010; Treu et al. 2010; Cardone & Tortora 2010) or the total 
mass profile (Cardone et al. 2009; Grillo 2010). On the other 
hand, one can also try to constrain the stellar IMF by mak- 
ing some assumptions on the dark matter content (Grillo & 
Gobat 2010). In the opposite regime, galaxy - galaxy lensing 
allows to investigate the outer regions estimating the total 
mass - to - light (M/L) ratio (Guzik & Seljak 2002; Hoekstra 
et al. 2005; Gavazzi et al. 2007) although large samples are 
needed to get reliable results. 

In spite of the quite large range spanned by their mor- 
phological and photometric properties, ETGs show several 
interesting correlations among their colors, luminosities, ve- 
locity dispersions, effective radii and surface brightness, the 
most famous example being the Fundamental Plane FP, a 
log -linear relation between the effective radius R e ff, the 
intensity I e = J(_R e //), and the central velocity dispersion 
ao (Djorgovski & Davis 1987, Dressier et al. 1987, Bender 
et al. 1992, Burstein et al. 1997). This fundamental plane 
is usually parameterized as: R e ff oc a^lt with (a,b) pre- 
dicted to be (2,-1) if ETGs are homologous systems in 
virial equilibrium and with a constant M/L ratio. The ob- 
served plane is, however, tilted with respect to the virial 
one since the different determinations of (a, b), depending 
on the photometric band and the sample used, are always 
different from the virial values. Jorgensen et al (1996) first 
derived a = 1.24 ± 0.07 and b = -0.82 ± 0.02 from a set 
of 225 early -type galaxies in nearby clusters observed in 
the r-band. While this result is consistent with the original 
observations of Djorgovsky & Davis (1987) and Dressier et 
al. (1987), it is nevertheless in striking contrast with the 
most recent determination relying on ~ 9000 ETGs ob- 
served within the framework of the Sloan Digital Sky Survey 
(SDSS). Using this large sample, Bernardi et al. (2003) have 
found a = 1.49 ± 0.05 and b = -0.75 ± 0.01, which are 
more similar to the K - band fundamental plane of Pahre et 
al. (1998). While the precise values of the FP coefficients 
are still debated, it is nevertheless clear that the observed 
FP is tilted with respect to the predicted plane. Such a 
tilt could be caused by a variation in the dynamical M/L 
ratio for ETGs as a result of a varying dark matter frac- 
tion (e.g., Padmanabhan et al. 2004; Boylan - Kolchin et al. 
2005, Cappellari et al. 2006, Tortora et al. 2009) or stellar 
population variations (e.g., Gerhard et al. 2001). Moreover, 
non - homology in the surface brightness profiles of elliptical 
galaxies (e.g., Graham & Colless 1997; Trujillo et al 2004) 
may be an other explanation of the FP tilt*. 



* Hereafter, with an abuse of terminology, we will refer to the 



All the above results have been obtained assuming that 
the classical Newtonian theory of gravity may be used also 
on galactic scales. However, the outer regions of galaxies 
typically are in a low acceleration regime and in this regime 
Newtonian dynamics has never been experimentally tested. 
Motivated by this consideration, Milgrom (1983) proposed 
to modify Newton's second law of dynamics as F = mg, 
where the acceleration g is now related to the Newtonian 
one gN as gfi(g/ao) — gN- Here, ao is a new universal con- 
stant and n(x) may be an arbitrary function with the prop- 
erties n(x » 1) = 1 and fi(x « 1) = x. The theory 
thus obtained, referred to as Modified Newtonian Dynamics 
(MOND), provides flat rotation curves for spiral galaxies 
and, as a by-product, gives a theoretical interpretation of 
the empirically determined Tully- Fisher law. Since its be- 
ginnings, MOND has been widely tested as an alternative 
to dark matter with remarkable success on galaxy scales 
(see, e.g., Sanders & McGaugh 2002, Milgrom 2008 and refs. 
therein) receiving a renewed interest after the proposal of a 
possible fully covariant relativistic formulation referred to as 
TeVeS (Bekenstein 2004). 

Motivated by these successful results, we model ETGs 
as a single component system derived from the observed 
luminosity profile and check whether the resulting MOND 
dynamics is consistent with the data. To this end, assum- 
ing a constant anisotropy model for the velocity dispersion 
profile, we compute the aperture velocity dispersion a ap 
and compare it with the observed value for a large sam- 
ple of local ETGs. Fitting the a ap data provides us a con- 
sistency check of MOND. Indeed, should the MOND pre- 
dicted values be smaller than the observed one or ask for 
unrealistic anisotropy profile, one can argue that DM is still 
needed. Should, on the contrary, MOND overestimates a ap , 
one should reconsider the choice of the interpolating function 
n(x) or the acceleration constant ao- As a byproduct, this 
test also allows us to estimate the anisotropy profile which 
can then be used to infer how the theoretical FP looks like 
in the MOND framework. It is worth noting, however, that 
a ap mainly probe the inner regions of the galaxy where dark 
matter is not expected to play a dominant role. As a conse- 
quence, although our analysis is definitely in the context of 
the MOND vs DM controversy, we do not expect to give a 
conclusive answer on this debate, but we nevertheless make 
a further step towards a possible solution. 

The plan of the paper is as follows. In Section 2, we 
present the mass density profile adopted for describing the 
ETGs luminous component, while the derivation of the aper- 
ture velocity dispersion in the MOND framework is given in 
Section 3. We then show, in Section 4, how MOND pre- 
dicts a FP-like relation for ETGs which we refer to as the 
MOND Fundamental Plane (MFP). Section 5 is then de- 
voted to both motivating the MFP and testing the viability 
of MOND in reproducing the observed dynamics of ETGs, 
while the predicted MFP coefficients are given in Section 
6 and compared to those of the observed FP. We finally 
summarize our results and discuss possible implications in 
Section 7, while in the Appendix A and B we briefly investi- 



plane with coefficients (a, b) = (2,-1) as the virial FP even if, 
as explained, deviations from this plane do not imply that ETGs 
are unvirialized systems. 
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gate the wavelength dependence of the MFP coefficients and 
discuss the impact of deviations from the main assumptions 
considered in the paper. 



2 THE PS MODEL 

In spite of their wide mass, size and luminosity range and 
the different chemical and stellar population characteristics, 
ETGs present a remarkable similarity in their photometric 
properties. As many studies show, (Caon et al. 1993; Gra- 
ham & Colless 1997; Prugniel & Simien 1997), their surface 
brightness is well described by the Sersic (1968) profile: 



I(R) = I e exp < — b n 



R 



R, 



l/n 



(1) 



with R the cylindrical radius^ on the plane of the sky and 
Ie the luminosity intensity at the effective radius Reff- The 
constant b n is determined by the condition that the lumi- 
nosity within R e f / is half the total luminosity. A very good 
analytical approximation is given by (Ciotti & Bertin 1999) : 



2n- 



0.009876 



n 



The deprojection of the intensity profile in Eq.(l) is straight- 
forward under the hypothesis of spherical symmetry, but, 
unfortunately, the result turns out to be somewhat involved 
combinations of the unusual Meijer functions (Mazure & 
Capelato 2002). In order not to deal with this cumbersome 
expression, we prefer to use the model proposed by Prugniel 
and Simien (1997, hereafter PS) whose three dimensional 
luminosity density reads : 



j(r) = jo 



R. 



eff 



with 



Jo = 



exp 



T(2n) 



-b n 



R, 



eff 



l/r, 



2R eff r[n(3-p n )] 



(2) 



(3) 



Here, F(a) is the F function, Iq = I{R — 0) = I e o hn , while 
the constant p n is chosen so that the projection of Eq.(2) 
matches a Sersic profile with the same values of (n, Reff, Ie)- 
A useful fitting formula is given as (Lima Neto et al. 1999; 
Marquez et al. 2001) : 



pn 



1.0 



0.6097 0.00563 
+ 



n 



In the following, we will be interested in the mass rather 
than luminosity density. Under the light - traces - mass as- 
sumption, the two quantities are immediately related as 



t Note that we have implicitly assumed that the intensity I does 
not depend on the angular coordinates. Actually, the isophotcs 
are not concentric circles, but rather ellipses with variable cllip- 
ticitics and position angles so that / = I(R, <p). However, in order 
to be consistent with our assumption of spherical symmetry of 
the three dimensional mass profile, we will neglect such an effect 
and, following a common practice, circularize the intensity profile 
considering circular isophothes with radii equal to the geometric 
mean of the major and minor axes. 



p(r) — T*j(r) with T, the stellar M/L ratio. For later appli- 
cations, it is convenient to define the following dimensionless 
quantity 

p(v) 

Peff 

with rj = r/R e ff and 

Mi, (,'• 



p(v) = 



V 



exp[-b n (r, 1/n -l)] , 



(4) 



Peff = P(V = 1) = T*j e~ 



47ri?3 // nT[n(3- Pn )] 



(5) 



Here M* = Y*Lt is the total stellar mass with Lt = 
27rnfe^ 2n e'' ,i r(2n)7ei?e// the total luminosity of the pro 
jected Sersic profile. Because of the assumed spherical sym 
metry, it is only a matter of algebra to show that the cumu 
lative mass profile for the PS model reads : 

7 [n(3-p„),M 1Al ] 



M(r) = 



r[n(3-p„) 



(6) 



where 7(0, x) is the incomplete 7 function. It is useful to 
introduce the following scaled mass profile : 



M{n) = 



Mjrj) = 7[n(3-p n ),M 1 H 
Meff 7[n(3-p„),&„] 

7[n(3 - p n ), bn 



with 

M eff = M{r) = 1) = Mi, 



r[n(3-p„)] 



(7) 



(8) 



As a final remark, let us stress that both M* and M e // 
may be determined from the measurement of the photomet- 
ric parameters (n, Reff, Ie) provided that an estimate of the 
stellar M/L ratio T* is available (for instance, from the re- 
lation between T* and the colours or from fitting the galaxy 
spectrum to stellar population synthesis models) . 



3 APERTURE VELOCITY DISPERSION 

A widely used probe to constrain the model parameters is 
represented by the line of sight velocity dispersion luminos- 
ity weighted within a circular aperture of radius R ap . This 
can be easily evaluated as : 



2 

^ an 



j" ap I(R)af os (R)RdR 
f* ap I(R)RdR 



(9) 



with ai os (R) the velocity dispersion projected along the line 
of sight. In order to compute this latter quantity, we first 
need to solve the Jeans equation for the radial velocity dis- 
persion a r . Following Sanders (2000) and assuming spherical 
symmetry, we write it as : 



dW)fr{r)\ + , 
ar r 



(10) 



with /3(r) = 1 — o'g/a r the anisotropy profile and g(r) 
the acceleration law. In the classical Newtonian dynamics, 
g(r) = GM t ot(r)/r 2 with M tot (r) the total (stellar + DM) 
mass. In the MOND framework, no DM is added so that 
M to t{r) = M(r) with M(r) the stellar mass from Eq.(6). 
Second (and most importantly), the acceleration g(r) is ob- 
tained by solving 



9(r)fJ. 



9(r) 

d 



: 9N{r) 



(11) 
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with fj,(x) the MOND interpolating function and ao = 1.2 x 
1CT 10 m/s 2 the MOND acceleration scale. Without loss of 



L 2 = 2nne D "b n M I e Riff^(2n, b n £l' p n ) , 



generality, we conveniently define : 
g(r) = iMOND{r)g N {r) 



(12) 



where the functional expression of jmond(v) will depend 
on the mass model and the /i(:r) function. Although fj,(x) = 
ir/Vl + x 2 has been the first proposal to be tested with suc- 
cess (Sanders & McGaugh 2002), recent analyses (Famaey 
& Binney 2005; Zhao & Famaey 2006; Famaey et al. 2007; 
Sanders & Noordermeer 2007; Angus et al. 2008b) seem to 
favour the simple form (Famaey & Binney 2005) 



H(x) = 



1 + x 



(13) 



which we will adopt in the following. Using Eq.(6) for the 
mass profile, inserting Eq.(13) into Eq.(ll) and solving for 
g(r), we finally find: 



2^MOND(ri) 

with 



1 + W1 + 



4q ?7 2 7[n(3 - p n ),b n ] 



a e//7["-(3~Pn),&n?? 1 / n ] 

GM eff _ GM* 7[n(3-p„),&„] 



R 2 efJ r[n(3-p n )] 



(14) 



(15) 



The computation of oi os may now be performed along the 
same steps as for the Newtonian case (see, e.g, Mamon & 
Lokas, 2005a,b) provided the following replacement rule 

GM(r)/r 2 -> GM{r)~fMOND(r)/r 2 

is applied. The final result turns out to be : 



I(R)o-L(R) 



2GM eff p eff /T* 



(16) 



with £ = R/R e ff, ij a — r a /R e ff a scaled anisotropy radius 
and K(r}/£,r]a,/ti) a kernel function depending on the choice 
of the anisotropy profile^. There are not many constraints 
on what the correct anisostropy profile can be. Moreover, 
the results in the literature have all be obtained under the 
assumption of Newtonian gravity and dark matter so that 
extrapolating them to the MOND framework is not moti- 
vated. As a first approximation, we will therefore consider 
models with constant anisotropy, i.e. /3(r) = /3, where we 
remember that j3 may range from — oo (for a system with 
fully tangential orbits) to 1 (for radial orbits only). For the 
anisotropic PS model in a MOND framework, we finally get : 

P)£d£, (17) 
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4irG Peff M eff R 2 eff 
T*L 2 7[ n (3 -Pn),b n ] J Q 



with 



poo 



^mond{ii) 



7[n(3 - 



l/r. 



T[n(3 - p n 



dr) , 



(18) 



t For some j3{r) profiles, the kernel is analytic and can be re- 
trieved from Appendix B of Mamon & Lokas (2005b). 



(19) 



and p denotes the set of parameters on which the integral 
depends. In the FP studies, the aperture velocity dispersion 
is usually referred to a circular aperture of radius R ap = 
Reff/8, i.e. £ ap = 1/8. Following common practice, we will 
set ao = <T ap ((,ap = 1/8) and note that, because of Eqs.(17) - 
(19), it is: 

2 



O"0 



s(n,I e ,Reff,T*,P) 



having noted that both p e / / and M e j / may be expressed as 
function of the photometric parameters (n, R e ff, le) and the 
stellar M/L ratio Y*. Investigating the shape of the function 
s(p) will be the aim of the next section. 



4 THE THEORETICAL MOND FP 

It is just a matter of algebra to show that 



2 GM„ & n e- i, "F(2n) s(p s ) 



Reff nr[n(3-p„)] 7 [2n,6„(l/8) 1 /»] 



with 



1/8 



(20) 



(21) 



s(ps) = / Zios{€,p s )£d4 

Jo 

and we have introduced the subset of parameters : 
p s = (n,a eff ,P) 

with a e ff is a function of (n, I e , R e ff, T*) through 
(M e ff,R e ff). It is computationally convenient to use 
a e ff/ao as parameter when estimating s(p s ) since this is 
the only way the full set of parameters enters the definition 
of the jMOND(rj) function in the integral. 

We can now make use of the definition of average effec- 
tive intensity to write : 

M* — T*Lt = 27rT*(/ e )i?e// 

where (I e ) and I e are related by (Graham & Driver 2005) 
</ e ) = n&- 2 V"r(2n)/ e . (22) 
Note that, using (I e ), the parameter o e // now reads : 
27rG(£ e )7[n(3-p n ),M 



(23) 



r[n(3-p„)] 

where we have defined the average mass surface density as : 

<S e ) = T,(J e ) . (24) 

Because of Eq.(24), Eq.(20) can be seen as a relation among 
(n, Reff, (S e )) which can be easily solved by introducing 
logarithmic units as : 

log Reff = 2 log a - log (E e ) - log s(n, (E e )) 

- log/(n)-log(27rG) (25) 

with 



f(n) 



6„e- i, ™r(2n) 



(26) 



nr[n(3 - p n )] 

Let us now assume that, over the parameter space actually 
covered by galaxies, the two terms s(n, (E e )) and f(n) may 
be well described as power law functions of their arguments. 
In such a case, Eq.(25) may be approximated as : 
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log R e f f = a log ao + b log (E e ) + c log n + d 



(27) 



with (a, b, c, d) parameters to be determined as described 
later. A guess for their values can be obtained by means of 
the following expected scaling relations : 

S(n,(E e )) ocn°*(E e ) 6s , f(n) oc n a " , 

so that one should find : 



a = 2.0 



b=-( 



+ 1) 



Actually, deviations of s(n, (E e )) and f(n) from exact power 
laws may be, in principle, compensated introducing a depen- 
dence on ao thus making a deviating from 2. We have there- 
fore left this parameter free in order to accomodate possible 
deviations from the canonical value induced by the change 
from the exact Eq.(25) to the approximated (27). 

As a final step, we now have to remember that both the 
stellar M/L ratio T* (Prugniel & Simien 1996; Cappellari et 
al. 2006; Tortora et al. 2009) and the Sersic index n (Caon et 
al. 1993; Graham & Guzman 2003; Mamon & Lokas 2005a) 
may actually be correlated with the galaxy total luminosity 
L (hereafter, we drop the label T for convenience). In a first 
good approximation, we can write : 



log T* = a* log L + , 
log n = a v log L + p v . 



(28) 
(29) 



Using then Eqs.(22) and (24) and putting together Eqs.(27), 
(28) and (29), we finally get : 

log R e f f = CtMFP log a + PmFP log (I e ) + 1MFP (30) 

with 



CtMFP 



Pmfp = 



yMFP = 



1 — 2(ba* + cov) 

(1 + a±)b + ca u 
1 — 2(ba* + cct v ) 

(ba* + ca„) log (27r) + + c/3„ + d 



(31) 



1 - 2(ba* + ca v ) 

which defines what we call the MOND Fundamental Plane. 

According to the above derivation, the MOND Funda- 
mental Plane§ (hereafter MFP) should be a perfect plane 
with no scatter. Actually, there are three main sources of 
scatter. First, Eq.(27) is an approximation for the exact 
Eq.(25). Thus, we expect that this works more or less well 
depending on each galaxy. In other words, the difference 
AR e ff between the solutions of Eqs.(25) and (27) will be a 
function of (n,R e ff, (E e )) and hence will change from one 
galaxy to another introducing a scatter around the MFP. As 
a second issue, one should take into account that the n - L 
and T* - L correlations are affected by an intrinsic scatter 
that propagates to the MFP. Finally, although not explic- 
itly denoted above, it is nevertheless clear from Eq.(18) that 
s(n, (E e )) also depends on the anisotropy parameter p. Since 

§ It is worth noting that the problem of the FP in MOND has 
been first discussed in Sanders (2000). However, in that paper, 
Sanders did work out a galaxy model looking for a solution of the 
Jeans equations for a polytropic system, while, here, we start from 
what we observe. See, also, Sanders & Land (2008) and Sanders 
(2010) for a similar approach. 



this quantity changes from one galaxy to another, this will 
induce deviations from a fiducial MFP defined by setting 
P to a reference value. The combination of all these effects 
make the galaxies scattering around the MFP thus giving it 
a non zero thickness. 



5 TESTING THE MOND SCENARIO 

The previous sections have demonstrated how to derive a 
fundamental plane in the MOND framework. However, the 
above derivation heavily relies on two main assumptions. 
First, we have assumed that the exact relation (25) may be 
replaced by the approximated expression (27). Second, wc 
are implicitly assuming that one is able to reproduce the 
dynamics of ETGs in the MOND framework without in- 
troducing any further dark matter component. In order to 
test this latter hypothesis against observational data, one 
can try to fit the observed aperture velocity dispersion in a 
large sample of local ETGs with well measured photometric 
properties. In such a case we can assume the three photomet- 
ric quantities (n, R e f f , to be known so that imposing 
that the observed and theoretically predicted a\ v are equal 
gives a constraint in the parameter space (T*,/3). We can 
moreover resort to stellar population synthesis codes to infer 
the stellar M/L ratio Y* from the galaxy colors thus finally 
ending up with an estimate of the anisotropy parameter p. 
Note that such a sample is also useful to test the first as- 
sumption and suggests what is the region of the parameter 
space (n, R e ff, (E e )) we have to explore to see whether our 
approximation holds for realistic ETGs systems. 



5.1 The data 

As a first step to carry on the approach detailed above, 
one has to assemble a sample of ETGs as large as possi- 
ble. To this aim, we have started from the NYU Value- 
Added Galaxy Catalog (hereafter, VAGC) which is a cross- 
matched collection of galaxy catalogs maintained for the 
study of galaxy formation and evolution (Blanton et al. 
2005) and mainly based on the SDSS data release 6 (Adel- 
man - McCarthy et al. 2007) . Among the vast amount of 
available data, we use the low - redshift (hereafter, lowZ) cat- 
alog of galaxies with estimated comoving distances in the 
range 10 < D < 150 h~ 1 Mpc. We refer the reader to Blan- 
ton et al. (2005) and the VACG website^ for details on the 
compilation of the catalog. Note that the lowZ catalog is 
actually updated only to the fourth ( Adelman - McCarthy 
et al. 2006) rather than the sixth SDSS data release thus 
covering an effective survey area of 6670 square degrees. 

From the lowZ catalog, we remove all the galaxies with 
unmeasured ao- This first cut leaves us with 43312 out 
of 49968 objects with magnitudes in the five SDSS filters 
u'g'r'i'z' . Based on the available data, we assemble an ETGs 
sample by imposing the following selection criteria : 

(i) 2.5 < n < 5.5 with n the Sersic index in the i' band 
and the upper end dictated by the code limit n = 6.0; 



http://cosmo.nyu.edu/blanton/vagc/ 
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(ii) R90/R50 > 2.6 (Shimasaku et al. 2001) with R f the 
Petrosian radii containing f% of the total luminosity; 

(iii) (Jo > 70 km/s since the measurement of the velocity 
dispersion for such small mass system could be unreliable as 
explained in Bernardi et al. (2005); 

(iv) (g' - r')- <g'-r' < (g' - r')+ with (<?' - r')± = 
pM r + q±8, M r the absolute magnitude in the r' filter and 
the parameters (p, q, S) have been tailored from Fig. 2 in 
Bernardi et al. (2005) where a different ETG sample has 
been extracted from the SDSS DR2 (Abazajian et al. 2004). 

The final sample thus obtained contains 9046 galaxies out of 
an initial catalog containing 49968 objects. It is worth noting 
that most of the rejected objects have been excluded by the 
first three cuts (retaining only 9105 entries), while the fourth 
cut only removes 59 further galaxies. This is reassuring since 
the last cut is somewhat qualitative and based on a different 
set of selection criteria (Bernardi et al. 2005). We then use 
the data reported in the lowZ catalog for the galaxies in 
the above sample to collect the quantities of interest. In 
particular, the average effective intensity (7 e ) is given by : 



in _ 6w dex[(M t -M s )/2.5] 
iU X 2-kRI 



(32) 



with dex(:r) = 10 x , Mt the galaxy absolute magnitude cor- 
rected for extinction, evolution and cosmological dimming, 
Mq the Sun absolute magnitude in the given filter 1 1 , while 
Reff is here expressed in kpc rather than arcsec. To this 
aim, we simply use : 

R eff (kpc) = R eff (arcsec) x D A (z)/206265 

with Da{z) the angular diameter distance in Mpc. 

The velocity dispersion reported in the lowZ catalog 
is measured within a fixed aperture Rsdss = 1-5 arcsec, 
while ao entering the MFP refers to an aperture of radius 
Rap = Reff/8. To correct for this offset, we follow j0rgensen 
et al. (1995, 1996) and set: 



obs lowZ w 

c>o = 0o x 



Rsdss 

Reff/8 



(33) 



with a ° w the value in the catalog and R e f / in arcsec here. 

5.2 Estimate of the stellar M/L ratio 

The lowZ catalog galaxies have been observed in five pho- 
tometric bands so that we can use the color information in 
order to infer their stellar M/L ratios (Tortora et al. 2009). 
To this aim, we start assembling a library of synthetic stel- 
lar population models obtained through the Galaxev code 
(Bruzual & Chariot 2003) varying the age of the population, 
its metallicity and time lag of the exponential star formation 
rate and assuming a Chabrier (2001) initial mass function 
(IMF). Then, we use the tabulated (it' , g' ,r' ,il , z') appar- 
ent magnitudes (corrected for extinction) of each ETG to 
fit the above library of spectra (suitably redshifted to ETG 
redshift) to the colours, thus getting the estimates of Y* for 
each galaxy in the catalog. Note that these values may be 



II We use Mq = (5.82,5.44,4.52,4.11,3.89) for the u'g'r'i'z' fil- 
ters respectively as evaluated from a detailed Sun model reported 
in www . ucolick . org/ ~cnaw/ sun . html 



easily scaled to a Salpeter (1955) or Kroupa (2001) IMF by 
multiplying by 1.8 or 1.125 respectively so that we can ex- 
plore other IMF choices**. We finally get an estimate of the 
total stellar mass M* simply as Y*Lt with the total lumi- 
nosity Lt taken from the catalog itself after correcting for 
the loss of flux due to the SDSS use of Petrosian magnitudes. 
Because of the errors on the colours, our T* estimates are 
affected by uncertainties difficult to evaluate. As a possible 
way out, one can use a Monte Carlo -like procedure gener- 
ating a set of colours from a Gaussian distribution centred 
on each mean colour and standard deviation equal to the 
colour uncertainty. Fitting the colours thus obtained to the 
synthetic spectra for each realization, one could generate a 
distribution of fitted parameters and take the median and 
median scatter of such a distribution as an estimate of Y* 
and its uncertainty (see Tortora et al. 2009 for further de- 
tails). We have tested that such a procedure gives errors on 
Y* of order 10%, but applying this method to the full sam- 
ple is too time demanding so that we prefer to neglect this 
source of uncertainty. Should we have used our data to con- 
strain the parameters of each single galaxy, this choice could 
have lead to an underestimation of the errors. However, we 
are here mainly interested in the statistical properties of the 
ensemble rather than fitting each individual galaxy. The un- 
certainty on each galaxy model parameters simply shifts the 
position of the galaxy in the parameters space, but does not 
alter the full distribution. Therefore, we are confident that 
neglecting the error on Y* has no impact on our results. 



5.3 Motivating the MFP 

Having estimated the stellar M/L ratio Y*, we are now able 
to derive for each galaxy in the catalogue the parameters 
(n, Reff, (S e )) we need to compute the central velocity dis- 
persion (To- We first evaluate the function s(p s ) over a finely 
spaced grid in (log n, log (a e ff/ao), j3). Considering typical 
values for the galaxies in the sample, we set the grid limits 
as: 2.0 < n < 6.0 and —2.5 < log(a e ///ao) < 2.5. Since we 
have no a priori information on what values the anisotropy 
parameter can take in a MOND scenario, we conservatively 
allow /3 to run in the range (0, 1) thus considering the full 
range for radial anisotropy. 

Having thus computed s(p s ), we can then resort 
to Eq.(20) to estimate cro over a finely spaced grid in 
(n, log Reff, log (S e )) for different values of /?. For each 
model in the grid, we then solve Eq.(25) and finally fit these 
values using Eq.(27) to estimate the coefficients (a,b,c,d) 
as function of the anisotropy parameter. In order to quan- 
tify the quality of the approximation, we also estimate 
A rms with A = log Reff (ex) — logR e ff(fit) and R e ff(ex), 
R e f f (fit) the exact and approximated solutions. 

It turns out that Eq.(27) indeed works quite well with 
0.06 < A rma < 0.09 and A values that are actually quite 
smaller than the rms one over most of the parameter space 
explored^. Moreover, there is no trend with any of the pa- 



** While this is correct for a Salpeter IMF, since it differs from a 
Chabrier IMF only for the low mass slope and predict very similar 
colours, for a Kroupa IMF this is not strictly true, but we could 
assume the scale factor above as a good approximation, 
tt As an alternative way of quantifying the goodness of the 
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rameters so that we can conclude that replacing the exact 
solution Eq.(25) with the approximated one Eq.(27) does 
not introduce any bias thus giving a theoretical support to 
our derivation of the MFP. 

Such a test makes also possible to infer the values 
of the (a, b, c, d) coefficients and how they depend on the 
anisotropy parameter. We first note that the use of the ap- 
proximated formula has made the coefficients (a, b) to de- 
viate from the values (2, —1) one should have inferred from 
Eq.(25) in the Newtonian case. While this is only a marginal 
difference for a, the effect is quite important for 6. Moreover, 
it is clear that, although trends of all the coefficients with 
the anisotropy parameter are clearly detected, the variation 
may be essentially neglected for a. On the other hand, both 
b and c significantly depend on P so that we expect a scatter 
in the MFP coefficients (umfp, Pmfp) due to the variation 
of the anisotropy parameter from one galaxy to another. We 
finally stress that the dependence of d on j3 does not have 
any impact on the slope of the MFP, but introduces a scat- 
ter in the zeropoint 'jmfp which translates into a scatter of 
galaxies around the best fit plane. 

One could wonder whether these results depend on the 
adopted MOND interpolating function. To this aim, we have 
repeated the above analysis considering the standard form 
for fj,(x) hence setting 

V \ a eff 7M3-P„),M 1/ n \J 

in Eq.(18). We find that the values of (a, b, c) are essentially 
the same thus suggesting a very weak dependence of the 
theoretical MFP coefficients on the particular fi(x) function 
considered. This can be qualitatively explained noting that 
(Jo is evaluated inside a very small aperture. Although for- 
mally the integral entering X; os (£) is defined along the full 
line of sight, integrating it over the very inner region of the 
galaxy, where 7mo at £>(?/) ~ 1, makes the details of this func- 
tion unimportant and explains why the choice of fJ,(x) has 
such a small impact on MFP coefficients. 



5.4 Matching the data 

While the above discussion shows that the derivation of 
the MFP is theoretically well founded, we still have to ver- 
ify that MOND is indeed able to match the dynamics of 
ETGs without resorting to dark matter. For each galaxy 
in the catalogue, we have both the photometric quantities 
(n,R e ff, {I e }) and the stellar M/L ratio T* so that, using 
Eq.(20), we can compute the aperture velocity dispersion as 
a function of the anisotropy parameter only. Matching the 
observed a ap to the theoretically predicted one makes then 
possible to constrain the constant anisotropy parameter p. 
Note that, to this end, we directly use the observed a ap , 
while <To is used when considering the MFP. Two consid- 
erations motivated this choice. First, we have estimated ao 
converting a ap through Eq.(33). Actually, this relation has 



approximation, one can note that the quantity A = 1 — 
log R e f f (fit)/ log R e f f (ex) is of order 0.01% over most of the 
parameter space explored. 



been derived assuming the validity of Newtonian gravity^ , 
while our analysis is in the MOND framework. Although it 
is likely that Eq.(33) still applies in the MOND case being 
do evaluated within the a/a » 1 region, it is a more con- 
servative choice not to use any correction and then checking 
a posteriori the validity of Eq.(33). As a second motivation, 
we note that ao is evaluated within 7? e ///8 where baryons 
are likely the dominant contribution and the Newtonian 
regime for acceleration holds. On the contrary, R ap /R e ff 
(with R ap = 1.5 arcsec) spans the 95% confidence range 
(0.16,1.58) with Ra P /Reff = 0.54 as median value. Kine- 
matical analysis in the Newtonian framework suggest that 
the dark matter mass fraction within ~ IReff is of order 
30% (for a Salpeter IMF). We therefore expect that a ap is 
more suited than ao to test the validity of MOND since it 
refers to a region where DM may play a not negligible role. 

A further remark is in order here before the fitting anal- 
ysis. First, both the photometric quantities and the stellar 
M/L ratio are known with an error and the same holds 
for the observed a ap . In order to take this into account, 
one could randomly generate a large ensemble of values for 
(n, R e ff, (Ie), T*, a ap ) according to a multinormal distribu- 
tion with central vector and covariance matrix set by the 
observed values. Solving for /3 each time, we can finally es- 
timate ft and its error from the distribution thus obtained. 
Because of the large ETGs sample we are dealing with, such 
a work is computationally quite expensive. However, we are 
here mainly interested in checking whether MOND is able 
to fit the data rather than deriving the exact value of the 
anisotropy parameter for every single galaxy. Moreover, we 
need the value of /3 or logr^a to estimate (qmfp, Pmfp) 
and the scatter in the MFP so that we are interested in its 
probability distribution function. For these reasons, we will 
neglect the errors on the observed quantities and simply take 
their central values for each galaxy in the sample. 

It is worth stressing that the shape of probability dis- 
tribution for p depends on the adopted filter, because both 
T* and the photometric quantities (n,R e ff, (I e )) are wave- 
length dependent. Moreover, the sample of galaxies we will 
fit is not the same for all the filters. In fact, for a given galaxy, 
the Sersic index n changes from one filter to another so that 
it is possible that a galaxy with n(i') > 2.5 has n(f) < 2.0 
in an another filter /, thus looking more similar to a disky 
system rather than an ETC one. Moreover, because of code 
failures or observational problems, some galaxies in the u' 
and z' filters may have n > 5.5 or an unmeasured absolute 
magnitude so that they have to be rejected from the sample. 
We stress, however, that more than 90% of the galaxies in 
the starting catalog are present in all the subsamples used to 
determine the P distribution in the different filters. Here, we 
discuss only the results for the fit in the i! filter thus using 
all the galaxies in the sample. The analysis in the remaining 
bands is similar. 



it Indeed, Eq.(33) can be derived by assuming a galaxy 
model and the anisotropy profile and then computing the ratio 
cr oos /cr Q owZ as function of the mass and anisotropy parameters. 
This is indeed how j0rgensen et al. (1995, 1996) derived their 
relation under the assumption that Newtonian gravity is correct. 
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5.4-1 Results for constant anisotropy models 

Since both the photometric parameters and the stellar M/L 
ratio have been fixed, we can straightforwardly solve the 
equation (n, R e f f , (S e ) , P) with respect to /3 

for each individual galaxy and then look at the distribution 
of the results. Such a procedure successfully works for all 
the galaxies in the sample with /3 values definitely pointing 
towards a strong radial anisotropy. This result is not sur- 
prising because, as well known in Newtonian gravity, radial 
anisotropy helps to lower the need for dark matter in the 
inner region of elliptical galaxies (Gerhard et al. 2001; De 
Lorenzi et al. 2008). For instance, De Lorenzi et al. (2008), 
using the NMAGIC code to fit the observed velocity disper- 
sion profile of NGC 4697, have found evidence of a mildly 
varying anisotropy profile with /3 ~ 0.3 in the inner regions. 
In the MOND framework, only part of the DM contribution 
is provided by the effect of the Jmond function entering a ap , 
while radial anisotropy supplies the remaining term needed 
to match the observed velocity dispersion. 

While radially anisotropic models are able to match the 
data, their physical consistency is not a priori guaranteed. 
In order to be self consistent, our density profile should be 
derived from a distribution function of the form L~ 2 ^ /e(E) 
with L the specific angular momentum and /e (E) a function 
of the binding energy E. The function Je{E) has to be pos- 
itive definite in order the model to be physically meaning- 
ful. Necessary conditions for the consistency of anisotropic 
models in the Newtonian framework have been derived in 
the literature (Ciotti & Pellegrini 1992; An & Evans 2006; 
Ciotti et al. 2009; Ciotti & Morganti 2010). Some of them 
may be extended also to MOND models taking into account 
that, for this purpose, MOND models behave like Newto- 
nian ones in an external effective potential ^ e // (Sollima 

6 Nipoti 2010). This latter can be derived by demanding 
that the MOND acceleration g(r) entering Eq.(10) equals 
d&eff/dr, i.e. one must integrate (numerically) the equa- 
tion d^eff/dr — Jmond(v)GM(< r)/r 2 with the condition 
that ^eff(r) vanishes at infinity. 

In the Newtonian framework, An & Evans (2006) have 
shown that, if the potential and the density are self con- 
sistently related, a necessary condition for the positivity of 
the DF of constant anisotropy models is 7 > 2/3, where 

7 = —d lnp/dlnr| r=0 is the logarithmic slope of the density 
profile at the centre. Ciotti & Morganti (2010) have however 
shown that this result can be extended to models with an 
external potential which is a case MOND can be compared 
to. We can therefore use the An & Evans result for the con- 
sistency of MOND constant anisotropy models. For the PS 
case, it is 7 = p(n) so that f3 < p(n)/2 must be taken as an 
upper limit on /3 to have physically consistent models. 

In order to take the above limit into account, we adjust 
both j3 and an effective M/L ratio kT*, where k accounts for 
deviations from the stellar M/L computed assuming a uni- 
versal Chabrier IMF. To this end, for each galaxy in the sam- 
ple, we determine k varying over the range (— 0.5,/3 max ) 
with p m ax = p(n)/2 and finally take as our best estimate 
the one with the smallest value of |1 — k\. Note that, in 
this way, we get models that are both physically consistent 
and with only minor deviations from the initial recipe for 
the stellar M/L ratio. Imposing 0.75 < n < 2.25 for the 
selected solution (see later for the motivation of such a se- 



lection criterium), we find that 92% of the galaxies may be 
successfully fitted with this procedure. Specifically, from the 
j3 and k distributions, we find : 

</3) = 0.07 , /3 med = 0.18 , 

68% CL : (-0.40,0.28) , 95% CL : (-0.45,0.37) , 

(k) = 0.98 , K med = 1.00 , 

68% CL : (0.96, 1.02) , 95% CL : (0.80, 1.03) . 

It is worth wondering whether the spread in k can be fully 
ascribed to variations in the stellar population properties. 
Indeed, the estimate of Y* from the galaxy colours depend 
on the details of the stellar population synthesis code we 
have used. There are many ingredients entering this code 
so that it is not fully unrealistic to expect that they can 
change from one galaxy to another, while here we have as- 
sumed them to be universal. For instance, a simple recipe 
to increase T* is to change the IMF from the Chabrier one, 
that we have used here, to the Salpeter one. It is worth not- 
ing that there is indeed still an open debate on what the 
IMF actually is. On one hand, direct star counts and obser- 
vations (mainly, in the Milky Way) point towards a Chabrier 
(or Kroupa) IMF, but our Galaxy is a spiral. On the other 
hand, gravitational lensing (see, e.g., Treu et al. 2010) and 
studies of the central DM fraction in local ETGs (Napolitano 
et al. 2010) suggest that also a Salpeter IMF can reconcile 
data and observations for ETGs provided a NFW model is 
assumed for the dark matter halo. Motivated by these con- 
trasting results, one can not exclude an IMF varying with 
luminosity (Renzini & Ciotti 1993; Tortora et al. 2009) or a 
correlation between the timescale of exponential star forma- 
tion rate and the luminosity. Moreover, the metallicity and 
the dust content may also change on a case -by -case basis. 
Considering a conservative ~ 25% uncertainty on T* for a 
given IMF and that T* has to be scaled by a factor 1.8 when 
replacing the Chabrier IMF with the Salpeter one, we finally 
get 0.75 < k < 2.25 as a range for this parameter. Should 
k be outside this (conservative) range, one should conclude 
that dark matter is needed to fill the gap between stellar 
and dynamical mass. Needless to say, the results we get for 
k are fully consistent with the statistical uncertainty on the 
estimated T*. While this outcome is expected by construc- 
tion, the fact that we are able to find a physically acceptable 
value of /3 makes us argue in favour of MOND being able to 
fit the velocity dispersion data with no dark matter provided 
a reasonable amount of anisotropy is allowed. 

It is worth noting that the confidence ranges for /3 are 
strongly asymmetric. Actually, 75% of the fitted models have 
f3 > with the remaining 25% giving rise to a long flat 
tail towards negative /3. If we simply cut our sample only 
retaining isotropic and radially anisotropic solutions, we find 
that the velocity dispersion may be fitted for 70% of the 
galaxies. Actually, this fraction is still larger if we repeat the 
above analysis, but now imposing /3 > 0. The fit turns out 
to be successful for 87% of the sample giving the following 
values characterizing the /3 and k distributions : 

</3)=0.17 , /3 med = 0.20 , 

68% CL : (0.10,0.30) , 95% CL : (0.0,0.37) , 
(k) = 0.98 , n med = 1.00 , 
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68% CL : (0.95,1.01) , 95% CL : (0.78,1.03) . 

The spread in j3 is now reduced with no statistical mean- 
ingful shift of the median value, while the shift in the mean 
is a consequence of the more symmetric distribution. The 
large fraction of successfully fitted galaxies then leads us to 
conclude that a mild radial anisotropy allows us to get phys- 
ically consistent models able to fit the data with no need for 
dark matter in a MOND framework. 

It is interesting to note that /3 does not correlate nei- 
ther with the photometric parameters (n, R e f / , (J e )) nor 
with the total luminosity L. Such a result has an important 
outcome. Since the MFP coefficients (omfp, Pmfp^mfp) 
given by Eq.(31) depend on /3 through (a,b,c,d), a strong 
correlation of /3 with L could generate a correlation between 
{ctMFP, Pmfp, Jmfp) and Lt and one should find different 
MFPs for different luminosity bins. However, this effect is 
likely to be weak because (a, 6, c, d) are actually weak func- 
tion of P; moreover, we need to have an MFP independent of 
luminosity because the observed FP is approximately unique 
for all the ETGs. 

As discussed in Appendix B, a correlation of k with L 
could change the MFP coefficients. On the one hand, such a 
correlation is expected noting that the dynamical M/L ra- 
tio Tdyn = kY* is found to correlate with luminosity when 
fitting data in the Newtonian + dark matter framework. Sur- 
prisingly, we do not find any scaling relations between k and 
L. Indeed, this is a consequence of how we select the best 
fit model for each galaxy. Let us suppose that a relation like 
k = k s (L/L s ) s exist for models with /3 fixed. Depending on 
(k s , L s , 8) and the galaxy luminosity, one could then obtain 
a value of k > 2.25. But we have chosen as best fit model 
the one having as small a value as possible of |1 — k\ so that, 
should k become too large, our algorithm would change /3 
(and hence k) to minimize 1 1 — k\ . As a consequence, unless 
the correlation is quite strong, our selection criterim washes 
out a possible scaling of k with L. We stress, however, that 
this is not a limitation of our analysis, but rather a conse- 
quence of consistently working in a MOND framework. 



5.4-2 MOND vs Newtonian gravity 

The above results convincingly show that MOND is able 
to fit the aperture velocity dispersion in ETGs without the 
need for additional DM. Actually, the data we are using 
mainly probe the region R < R e ff which is not in the deep 
MOND regime (a << ao). Indeed, for the galaxies in the 
sample, the median value of log(a e ///a ) is 0.56, while the 
68% confidence range is given by 0.21 < log(a e ///ao) < 
1.87, i.e. 1.6 < a e ff/ao < 7.4. We are therefore in an in- 
termediate regime where the total acceleration a is roughly 
of the same order of magnitude as ao. We therefore expect 
that it is possible to fit the same data without DM also in 
the Newtonian regime. To this end, we have used the same 
equations as above setting 7MOArr>(r) = 1 to recover the 
Newtonian formulae. Moreover, we have still to take into 
account the upper limit on /3 so that we use the same pro- 
cedure as above adjusting the anisotropy parameter and the 
ratio k between the dynamical and the stellar M/L ratio. 

Imposing the same cut on k, we find that ~ 80% of the 
galaxies can be successfully fitted, a fraction smaller than in 
the MOND case, but still satisfactorily large. It is, however, 



interesting to note that, for most of the galaxies, the value 
of j3 is quite close to the upper limit, while typical n values 
are larger than in the MOND case. Indeed, we get : 

</3)=0.33 , /3 med = 0.40 , 

68% CL : (0.02, 0.42) , 95% CL : (0.0, 0.44) , 

(k) = 1.35 , K med = 1.29 , 

68% CL : (1.00,1.79) , 95% CL : (0.81,2.13) . 

As it is apparent from the k confidence ranges, one must in- 
crease the M/L ratio up to very large values that, while still 
marginally consistent with deviations from the stellar one 
because of a different IMF, may be also easily interpreted in 
terms of a non negligible dark matter content. Needless to 
say, this is consistent with expectations since we know that, 
unless an unreasonably high radial anisotropy is introduced, 
Newtonian gravity can not explain ETGs dynamics without 
the boost provided by a dark matter halo. 

It is worth comparing the values of (/3, k) obtained in 
the MOND and Newtonian cases for each galaxy. We find : 

((1 + /3m)/(1 + /3jv)> =0.85 , /3 med = 0.86 , 

68% CL : (0.76, 0.95) , 95% CL : (0.71, 1.00) , 

(km/kjv) = 0.64 , K med — 0.67 , 

68% CL : (0.41,0.85) , 95% CL : (0.24,0.93) . 

These values clearly show that Newtonian models have a 
larger radial anisotropy and k value than MOND ones in 
agreement with the expecations. Indeed, for given model 
parameters, the boost to the Newtonian velocity dispersion 
provided by the ^mond^) factor reduces the need for both 
radial anisotropy and deviations from the fiducial stellar 
M/L. Not surprisingly, we find that km/kn is correlated 
with log(a e ///oo) so that the smaller is this quantity, the 
smaller is km/kn- This is exactly what one can anticipate 
noting that the smaller is log (a e ///ao), the larger is the re- 
gion of the galaxy where the MOND regime applies so that 
the smaller is the value of km needed to fit the data for the 
same value of the anisotropy parameter. 



6 THE MFP COEFFICIENTS 

Eq.(31) makes it possible to estimate the MFP parameters^ 
(ctMFP, Pmfp) provided the values of (a*, a„) are given and 
the anisotropy profile parameter has been set so that (a, 6, c) 
may be computed. The analysis in the previous section has 
demonstrated that MOND is able to fit the measured aper- 
ture velocity dispersion of the assembled ETGs sample and 
has enabled us to determine the distribution of /3 for the 
galaxies in our sample. In order to compute the MFP coef- 
ficients, we first need to determine (a*, a„) as we discuss in 
the following subsection. 



3S Hereafter, we will parameterize the MFP through the two 
slope related quantities (cimfp, 0m fp) without considering any- 
more the zeropoint -fMFP- 
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Table 1. Calibration parameters and intrinsic scatter for the luminosity scaling correlations in Eqs.(28) and (29) between logL and 
logY* and logn respectively. Columns are as follows: 1. filter id; 2. maximum likelihood parameters (q*, /?*, <Tj n t) f° r the logL-logY* 
relation; 3, 4. median value and 68 and 95% confidence ranges for (a*,<7; n t); 5. maximum likelihood parameters (a v ,fSv, (Tint) for the 
log L- logn relation; 6., 7. median value and 68 and 95% confidence ranges for (a„,<7j n t). 



Id 



(a*, 0*,<Tint)ml 



(a*)tl 



+1ct +2a 
2a 



la + 2a 
2a 



{ot v ,P v , (Tint) ml 



+ la +2a 
2a 



((Tint)t\ 



la +2a 
2a 



(0.0,0.680,0.0) 0.0 

(0.037,-0.115,0.0) 0.037 
(0.021,0.115,0.0) 0.022 
(0.004,0.290,0.0) 0.005 

(-0.012,0.406,0.0) -0.010 



+0.202 +0.416 
-0.201 -0.416 



+0.176 +0.276 
-0.135 -0.235 



+0.129 +0.264 
-0.128 -0.262 



+0.112 +0.228 
0.112 -0.228 

+0.125 +0.265 
-0.145 -0.285 



0.060 



+0.077 +0.211 
-0.043 -0.057 



fl nzW+°- 054 +0 144 
u.u<±o_ 030 _ 041 

n n , -,+0.053 +0.140 
u - u ^ -0.029 -0.039 



0.039 



+0.048 +0.125 
-0.028 -0.037 



n nAq+ - 054 +0-143 

u.u<iO_ Q31 _ .04i 



(0.039,0.146,0.0) 0.03!) 
(0.050,0.025,0.0) 0.050 
(0.047,0.075,0.0) 
(0.052,0.047,0.0) 
(0.044,0.114,0.0) 0.045 



+0.130 +0.270 
-0.131 -0.192 



+0.139 +0.237 
-0.093 -0.136 



n rvL7+°- 088 +0-180 

u.u<il _ _o88 _o.l79 

n n^9+ 078 +0-159 

u - udz -0.077 -0.157 



+0.094 +0.192 
-0.094 -0.191 



0.041 



+0.054 +0.145 
-0.029 -0.039 



n rm +0.040 +0.106 
u.uoi_ 022 _ 031 



0.030 



+0.047 +0.099 
-0.021 -0.047 



fl n9Q+°- 055 +0-092 
u.uz»_ Q 021 _ 055 



0.031 



+0.039 +0.103 
-0.022 -0.031 



6.1 Luminosity scaling relations 

Eqs.(28) and (29) are linear relations (although in a loga- 
rithmic space) so that to determine their parameters we can 
resort to a general Bayesian procedure described in detail in 
D'Agostini (2005). Let us suppose that (R, Q) are two quan- 
tities related by a power - law relation as R = BQ A with a 
certain intrinsic scatter ai„t ■ In logarithmic units, this reads 
log R = ai log Q + /3; with ai — A and /3i = log B. In order 
to determine the parameters (a, 13, (Tint), we then maximize 
the following likelihood function C = exp [— L(ai, fit, (Tint)] 
with : 



L(ai, Pi, (Tint) 



2 a \ 



2 



a t Xi 



a: 



2 

int 



(35) 

'Vi T "( u Xi 

where (xi,yi) = (log Qi, log Ri) and the sum is over the AT 
objects in the sample. Note that, actually, this maximization 
is performed in the two - parameter space (ai,(Ti n t) since Pi 
may be estimated analytically as : 

-l 

Vi 



UlXi 



£7 



so that we will no longer consider it as a fit parameter. The 
median values and confidence intervals for a given quantity 
can then be determined by studying the shape of the corre- 
sponding marginalized likelihood, i.e. the integral of £ over 
the other parameter. 

Let us then determine a* and a v using the method out- 
lined above. Note that these quantities depend on the wave- 
length adopted, but not on the anisotropy profile or MOND 
interpolating function. As a preliminary task, we bin the 
galaxies in roughly equally populated luminosity bins and 
remove all the galaxies with n(f) £ (2.0, 5.5) and problems 
with the absolute magnitude estimate. For a given luminos- 
ity bin and filter, we then analyse the distributions of Y* 
and n and assign to that bin a value for T* and n using 
the median as central value and the 68% confidence range 
as uncertainty. Note that, since such distributions may also 
be asymmetric, the errors on Y* and n may turn out to be 



asymmetric. However, Eqs.(35) and (6.1) assume that the 
uncertainties are symmetric. We therefore follow D'Agostini 
(2004) and correct the observed values as : 



Vc 



yobs + (A+ — A_ 



(A++A_)/2 



with y b s the median value, {y m in,ymax) its 68% confidence 
range and A + = y max — y a bs, A- = y obs — y m in- It is worth 
stressing that such a correction is actually quite small so that 
we are confident that it is not biasing anyway the estimate 
of (a*, a K ). The fit results are summarized in Table 1, where 
we report the values of the calibration parameters and the 
intrinsic scatter for the different filters. 

As a general result, we find that both a* and a„ are 
quite small and, within the large error bars, are compati- 
ble with zero. Even if we consider only the best fit values, 
it is nevertheless clear that the correlation is quite shallow. 
However, the wide confidence ranges on (a*,a„) prevent us 
from drawing any definitive conclusion other than the quali- 
tative observation that T* and n only weakly correlate with 
the luminosity. While the result for the Y* - L correlation 
is in agreement with previous results in the literature, the 
very shallow slope of the n - L relation is somewhat at odds 
with previous findings. For instance, using galaxies with well 
determined photometric parameters in the Virgo sample, 
Nipoti et al. (2008) found n oc L% 27±(h02 which is defini- 
tively larger than our best fit a v for the g' filter (the closest 
to the B one), although well within the wide 68% confi- 
dence range. Such a discrepancy is actually illusory and is 
due to the combination of the small luminosity range probed 
and the large error bars. Indeed, our ETGs sample actually 
spans approximately just one order of magnitude in luminos- 
ity (the first bin being at 9.1 and the last at 10.4). Adopting 
the Nipoti et al. slope, this leads to a change in n of order 
0.35 which is lower than the width of the 68% confidence 
range for n in a given bin. Indeed, if we fix the slope of the 
n - L relation in the g' filter to the one given by Nipoti et 
al. and only adjust the zeropoint and the intrinsic scatter, 
we get a very satisfactory result. We therefore argue that no 
discrepancy is actually present. As an alternative, we could 
have fitted the n-L relation without dividing galaxies into 
luminosity bins. However, binning allows to wash out the 
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Figure 1. Histogram of the hmfp values obtained in the V filter. 

biases in the n determination which are inherited by the 
SDSS photometric fitting code so that we have preferred to 
use this procedure even at the price of having larger error 
bars. 

6.2 Predicting the MFP coefficients 

The (cxmfp, /3mfp) values may be estimated using Eq.(31) 
provided one has set the filter (and hence the corresponding 
value of a* and a v ) and evaluated the (a,b,c) coefficients. 
However, we do not have a single value for all the quan- 
tities of interest, but rather a distribution that, as a first 
order approximation, we can model as a Gaussian one with 
central value and standard deviation obtained by correct- 
ing the asymmetries in the inferred distributions. In order 
to estimate ((Xmfp, Pmfp), for a given combination of fil- 
ter, anisotropy profile and MOND interpolating function, 
we generate (/}, a*, a v ) from a multinormal distribution and 
use Eq.(31) to compute the MFP parameters. We repeat 
this procedure 10000 times and use the distribution thus 
obtained to estimate the mean and median values and the 
68 and 95% confidence ranges. To this end, we use the p 
distribution estimated as described in Sect. 5.4.3 since that 
procedure allows us to both get physical consistent models 
and minimize the deviations from the fiducial T* value. 

The resulting cxmfp and Pmfp distributions for the i' 
filter are shown in Figs. 1 and 2 and can be quantitatively 
summarized by the mean and median values with confidence 
ranges given below : 

{(Xmfp) = 1-62 , {(XMFp)med = 1-58 , 

68% CL : (1.31,1.94) , 95% CL : (1.13,2.42) , 

{Pmfp} = -0.70 , (p M Fp)med = -0.69 , 

68% CL : (-0.74, -0.66) , 95% CL : (-0.80, -0.63) . 

As it is apparent from Figs. 1 and 2, while the omfp dis- 
tribution is quite wide, the Pmfp one is, on the contrary, 
quite narrow thus leading to small errors on this quantity. 
Such a result can be explained by noting that, since (a*, a v ) 
are quite small, Eq.(31) approximately reduce to cxmfp — ci 
and Pmfp — b. Since a depends more strongly on P than on 
b, the scatter of P around the median value has a stronger 
impact on cxmfp than on Pmfp thus explaining why the dis- 
tribution of cxmfp is wider than the distribution of Pmfp- 
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Figure 2. Same as Fig. 1 but for the Pmfp parameter. 

Both the mean values and the 68% confidence ranges 
are almost unchanged if one replaces the simple MOND in- 
terpolating function with the standard one. This result is 
encouraging since it tells us that the arbitrariness in the 
choice of the simple or standard function does not bias the 
final MFP coefficients. Actually, such a conclusion should 
be better tested by considering other /j,{g/ao) expressions 
leading to different 7mojv £>()") to be inserted in Eq.(17). We 
can, however, argue that the effect will be quite weak as can 
be easily explained by looking at the two "/mond{v) func- 
tions in Eqs.(14) and (34). Both expressions may be roughly 
approximated by jmond — 1 + (ao/fle// )/(??) with the de- 
tails of the function f(rf) depending on the adopted expres- 
sion of fj,(g/ao). For typical values of the ETGs parameters, 
do/cteff « 1 so that the details of the f{rj) function are 
important only for very large values of r\. However, the con- 
tribution to (Jo mainly comes from the inner regions thus ex- 
plaining why the choice of the MOND interpolating function 
only negligibly affects the estimate of the MFP coefficients. 

In the Newtonian framework, an easy application of 
the virial theorem (and the hypotheses of homology and 
constant M/L ratio) leads to what we have referred to as 
the virial plane, namely a FP - like relation with coefficients 
(a V ir,p v ir) = (2,-1). The MFP coefficients are definitely 
different from the virial ones. This is only partly due to the 
use of the scaling relations (28) and (29). If we turn off their 
effect by forcing a* = ct v — 0, we should have obtained 
values for {cxmfp, Pmfp) quite similar to those reported 
here. It is actually the combined effect of the anisotropy 
and MOND that leads to this departure from the Newto- 
nian virial predictions. The effect is however quite weak on 
cxmfp which only mildly departs from a V i r = 2 remain- 
ing consistent with this value within the errors, while it is 
stronger for Pmfp with the virial value p v i r = — 1 definitely 
out of the 95% confidence range. We therefore argue that 
MOND is indeed able to introduce departure from the New- 
tonian virial plane even if the ETGs inner regions are for 
the most part in the Newtonian regime. 

The disagreement between the virial FP coefficients and 
the observed ones is usually referred to as the problem of the 
FP tilt. The introduction of luminosity scaling relations for 
n and Y* is unable to explain the tilt so that one postulates 
the presence of a dark matter halo representing ~ 30—50% of 
the mass within R e ff and providing a total M/L ratio scal- 
ing with luminosity as T to t oc L° T with qt — 0.15 — 0.25. 
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It is worth wondering whether the FP tilt may instead be 
explained by our MOND based models. Should this be the 
case, the MFP coefficients should turn out to be consistent 
with the observed FP coefficients. Actually, there are dif- 
ferent estimates of the FP coefficients relying on different 
samples and different fitting procedures so that a fair com- 
parison is not easy (see, e.g., Bernardi et al. (2003) for a table 
with some estimates). We have therefore used our Bayesian 
method (generalizing from 2 to 3 parameters the previous 
formulae) to find as best fit parameters : 

{a obs , p obs , a int ) = (1.13,-0.75,0.081) 

in good agreement with (a obs , /3 obs , o int ) = (1.14 ± 
0.04, -0.76 ±0.01, 0.085) in Bernardi et al. (2003) which use 
a direct fit method and correct for evolution and selection 
effects. The median value and the 68% and 95% confidence 
ranges read^ : 



dobs 



n n R +0.33 +0.69 
u - yD -0.30 -0.52 



Pob 



n R ,+0.14 +0.25 
-U.DO_Q.21 -0.34 



The comparison of the predicted MFP coefficients with the 
observed ones shows that the tilt of the FP has been only 
partially reduced. Indeed, if we consider the best fit values 
for (a obs , /3 obs ), we find that, while /3 obs is nicely close to the 
median Pmfp, a obs is definitely smaller than the median 
Q.MFP- However, in a Bayesian framework, what is more im- 
portant is the comparison with the marginalized confidence 
ranges. In this case, Pm fp and /3 b s are in remarkable good 
agreement, while the 95% confidence ranges of cimfp and 
a b s have only a small overlap. We therefore conclude that 
there is only a mild chance that a MOND based approach 
align the observed FP with a theoretically predicted one. 

An important remark is in order here. We have shown, 
in Sect. 5.4, that MOND is consistent with the observed 
ETGs dynamics, that is to say, given (n, R e ff , (I e )) and 
the stellar M/L ratio T*, we can find a value for /3 such 
that the theoretically a ap equals the observed one. We have 
checked that this condition is enough to guarantee that the 
same matching still works for ao with ao estimated from a ap 
through Eq.(33). Since a galaxy is placed on the observed 
FP according to its a value, predicting the correct ao en- 
ables us to correctly place the galaxy on the observed FP. 
Nevertheless, our theoretical MFP is tilted with respect to 
the observed one. This inconsistency originates from the fol- 
lowing reason. When one considers an individual galaxy, the 
value of p is set according to the values of (n, Reff, T*) 
of that particular galaxy. On the contrary, the theoretical 
MFP is derived from a model that is assumed to be the 
same for all the galaxies. In particular, the value of /3 used 
to get (omfp, Pmfp) is the same for all the galaxies with 
the error on the MFP coefficients deriving from the width 
of the P distribution. Of course, P is not the same for all 
the galaxies so that one can not use the theoretical MFP to 
predict the ao value for a given galaxy. In other words, one 
could always find (omfp, Pmfp) for each individual galaxy; 
however, the set of these values does not define a single MFP, 



«^ Note that, to save computer time, we have fitted not the full 
sample, but a subsample made out by 1000 randomly selected 
ETGs. While this has no effect on the best fit parameters (as we 
have explicitly checked), it is likely that the 95% confidence range 
is wider than what we would have obtained using the full sample. 



but rather a set of N MFPs different from each other. The 
fact that this collection of planes does not reduce to a single 
one is another way of saying that the observed FP and the 
theoretical MFP are tilted. 



7 CONCLUSIONS 

The MOND framework has nowadays a long history of suc- 
cesses when applied to galaxy scale systems: it was shown 
to be able to efficiently explain spiral galaxy flat rotation 
curves, and empirical scaling relations such as the Tully- 
Fisher law and the Baryonic Tully - Fisher relation (see 
Sanders & McGaugh 2002 and refs. therein). In the stan- 
dard Newtonian scenario, dark matter is invoked not only 
in spiral galaxies, but also in ETGs to both reproduce the 
observed dynamics and explain the FP tilt. Since MOND 
is, by construction, a universal theory, one should be able 
to remove dark matter in ETGs as in spiral galaxies. Mo- 
tivated by this consideration, we have investigated whether 
this is the case and derived an FP - like relation which we 
have referred to as the MOND fundamental plane (MFP). 

To this aim, we have extracted from the NYU-VAGC 
lowZ catalog a sample of ~ 9000 ETGs with accurate five 
bands photometry and measured aperture velocity disper- 
sion. Assuming that no dark matter is present, we can sim- 
ply model an elliptical galaxy using the Prugniel and Simien 
(1987) model, setting the (n, Reff, {le}) parameters from 
photometry and inferring the stellar M/L ratio T* from 
the observed colours. We then consider a constant velocity 
anisotropy profile and the simple form for the MOND in- 
terpolating function /j,(g/ao) and constrain the anisotropy 
parameter f3 by matching the predicted and the observed 
aperture velocity dispersions. As a first important result, 
we find that physically consistent MOND models can be 
found to fit the data provided a radial velocity anisotropy 
is assumed and the dynamical M/L ratio is adjusted within 
the uncertainties of the stellar M/L. These results are in- 
dependent of the choice of the (poorly constrained) MOND 
interpolating function n(x) so that we may safely argue that 
MOND can remove dark matter not only in spiral galaxies, 
but also in early -type ones. 

Fitting the velocity dispersion data also allows us to in- 
fer the distribution of the anisotropy parameter and hence 
estimate the MFP slope coefficients (cumfp, Pmfp)- We find 
that these quantities are different from those expected when 
applying the virial theorem (with homology and constant 
M/L ratio assumptions) in the Newtonian framework, i.e. 
MOND is able to tilt the virial plane without the need of 
assuming a varying dark matter content. Such a tilt is, how- 
ever, still too small to align the MFP with the observed 
FP. Indeed, we find that, while Pmfp is in remarkable good 
agreement with Pots, cimfp is significantly larger than a b s 
so that we need a mechanism to further tilt the MFP. 

In particular, one could explore deviations from the 
pure MOND scenario we have used throughout the paper. 
Two possibilities are briefly hinted at here and discussed 
in some detail in Appendix B. As a first scenario, a vary- 
ing M/L ratio can reconcile the predicted MFP with the 
observed FP provided that K e ff = Y e ///Y* oc L 0,26 . A sim- 
ilar scaling might be expected in the MOND model by An- 
gus (2009) with 11 eV sterile neutrinos. Angus et al. (2010) 
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show that, in this model, more massive galaxy clusters, that 
generally harbour the most massive and luminous central 
galaxies, have larger central densities of sterile neutrinos 
and, presumably, a larger total M/L ratio. 

As a totally different approach, one can rely on the ex- 
ternal field effect (EFE). Different from the varying M/L 
scenario, the EFE is naturally motivated because it is a 
consequence of the nonlocal feature of the MOND theory. 
The inclusion of EFE is able to tilt the MFP without the 
need of additional matter and without resorting to some fine 
tuned mechanism. The scaling of g ext /oo with L needed to 
reconcile the theoretical MFP with the observed FP is rea- 
sonable, but has to be verified by a careful analysis. To this 
end, an ideal approach could be fitting the full velocity dis- 
persion profile ai os (R) rather than the aperture value a ap - 
This test has to be carried on for a statistically meaningful 
sample of ETGs with ai os (R) measured both in the inner 
regions (where the EFE has a small impact) and in the out- 
skirts (where the anisotropy profile reduces to its asymptotic 
value) so that the fi-g e xt degeneracy is broken. Moreover, 
this sample should cover a wide range in luminosity and 
probe different environments. 

As a final comment, we would like to stress that, al- 
though the two scenarios proposed above are, at the mo- 
ment, only speculative, reconciling the theoretical MFP with 
the observed FP and reproducing other ETGs scaling rela- 
tions in the MOND framework can open a new and fruitful 
way towards discriminating between modified dynamics (or, 
more generally, modified gravity theories) and dark matter. 
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APPENDIX A: WAVELENGTH DEPENDENCE 

Eq.(3f) shows that the MFP coefficients depend on the 
slopes (a*, a„) of the T* - Lt and n- Lt relations which are 
different depending on the filter considered, as reported in 
Table f. Moreover, as explained in Sect. 5.4, the distribu- 
tion of the anisotropy parameter also depends on the fil- 
ter adopted because (n,R e ff, (J e ),T*) vary with the filter. 
As a consequence of these effects, the MFP coefficients will 
depend on wavelength. To check how they change with A, 
we report in Table Af the values of (<xmfp,Pmfp) for the 
fiducial case of constant anisotropy with the simple MOND 
interpolating function. 

ft is immediate to see that the MFP coefficients are 
essentially the same (within the errors) whatever is the 
adopted filter. This is an expected result since both (a*, a„) 
and the distribution are very weak function of wavelength. 
This is in net contradiction with the observed FP coefficients 



Table Al. The MFP coefficients for the fiducial case in different 
filters. Columns are as in Table 2, while the rows are for filters 
«' , g' , r', i' , z' respectively. 
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1,ou -0.31 -0.52 


-0.71 


n fin+°- 04 + ' 07 

u - uu -0.06 -0.14 



that move towards the virial one (2, —1) as the wavelength 
increases. As a consequence, the tilt of the MFP changes 
with the filter thus leading further constraints on the un- 
derlying unknown phenomenon that has to be invoked to 
reconcile theory with observations. 



APPENDIX B: DEVIATIONS FROM A PURE 
MOND SCENARIO 

Up to now, we have adopted a fully MOND framework. 
Specifically, we have assumed that the visible matter is the 
only source of gravity generating the motion of the stars and 
that their acceleration may be computed from the Newto- 
nian one with Eq.(ll). As we have seen, this model is able to 
fit the ETG velocity dispersion data, but leads to an MFP 
which is still tilted with respect to the observed one. We 
therefore qualitatively discuss here what could be the im- 
pact of relaxing one of the two assumptions above to see 
whether deviations from our idealized MOND scenario may 
help to reconcile the theoretical MFP with the observed FP. 
We admit that this section is quite speculative, but we in- 
clude it here in order to propose some scenarios that could 
be tested with future data or a different tracer of the ETGs 
gravitational potential. 

Bl Varying M/L ratio 

According to the original idea motivating the introduction of 
MOND, nothing but the visible matter should be considered 
when modelling a galaxy and trying to fit the data. As a 
corollary to this assumption, the total M/L ratio T to t must 
equal the stellar T*. Even if in order to avoid too much 
radial anisotropy we have allowed to rescale the M/L by 
the parameter k, we have finally chosen the solution with 
k as close as possible to 1 to be consistent with the idea of 
no dark matter. Some general considerations, however, may 
be invoked to depart from the T to t — T* ansatz. First, it 
is well known that MOND needs a dark matter component 
in order to be in agreement with data on the galaxy cluster 
scale (Aguirre et al. 2001; Sanders 2003; Pointconteau & Silk 
2005; Angus et al. 2007, 2008a). Should this further term 
be represented by a 2eV neutrinos (Sanders 2003, 2007) 
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or 11 eV sterile neutrinos (Angus 2009; Angus et al. 2010) 
arranged in a hot dark halo, the total M/L ratio should be 
larger than the stellar one. On the other hand, our estimated 
T* depends on the ingredients used to obtain the starting 
library of stellar population models. Changing one of these 
ingredients would alter Y* possibly leading to an effective 
Ttot larger than our assumed stellar M/L value. As a simple 
way to account for this possibility, we redefine the effective 
surface mass density as : 



<E e ) = T tot (7 e ) = /s B T*(/ e > 



(Bl) 



If we assume that ks correlates with luminosity and that 
the log ke - log L may be well approximated by : 



log ke = & K log L + /3 K , 



(B2) 



we can easily rederive the MFP to show that Eq.(30) still 
holds provided the coefficients are redefined as : 



OiMFP = 



1 — 2arf> — 2a„c 
(1 + ar)b + ct v c 




1 — 2qt& — 2a„c 

(arb + a v c) log (2tt) + /3y + b + /3„c + d 
1 — 2arb — 2a v c 



(B3) 



(B4) 



Needless to say, introducing ke changes the distribution 
of the fitted anisotropy parameter and the values of the 
(a, b, c, d) coefficients. Assuming, as a first approximation, 
that the change in /3 is not too large, we solve olmfp = cx-obs 
with respect to a K . By using the median values of (/?, a t ,a„) 
for the constant anisotropy model and a obs — 1.13, we find 
a K ~ 0.35. It is worth noting that a similar scaling is ob- 
tained when comparing stellar and dynamical masses in a 
Newtonian stellar+dark halo framework (Padmanabhan et 
al. 2004; Tortora et al. 2009; Cardone et al. 2009). Actu- 
ally, reproducing such a scaling of kj with luminosity can 
be problematic in a MOND scenario. As quoted above, in 
order to have ks / 1, one has to postulate the presence of 
a halo made of (2eV or lleV) neutrinos. We now find that 
such a term should provide a contribution to the velocity 
dispersion that increases with the total luminosity in a sim- 
ilar way as the dark halo in Newtonian gravity. Explaining 
how neutrinos in the halo interact with the baryons in the 
stellar component in such a way to reproduce the needed 
scaling of ks with L can be a difficult hurdle to overcome. 

Should such a mechanism be found, the problem of the 
tilt of the MFP is indeed solved. Introducing a K — 0.35 
in the expression for Pmfp and using the median values 
of (/3,a*,a„) gives Pmfp — —0.64 which is essentially the 
same as for the kt. = 1 case. This can be explained by look- 
ing at the approximated expression reported above showing 
that the same term arb enters both the numerator and the 
denominator thus weakening the correction. As a result, the 
predicted Pmfp is still in agreement with the observed one 
thus completely solving the tilt problem. 



B2 The external field effect 

As a second possibility, we consider modifying the relation 
between Newtonian and MOND acceleration. We note that 
Eq.(ll) implicitly assumes that the galaxy is an isolated 
system. Since MOND is a non-local theory, the motion of 
the stars in the galaxy is actually determined not only by 
the galaxy potential, but also by the distribution of mat- 
ter outside the galaxy itself. This is usually referred to as 
the external field effect (EFE) and is taken into account by 
replacing Eq.(ll) with the following one: 



f 9+ 9ext 

V a 



9 = 9n 



(B5) 



where g ext is a constant (with typical values ~ 0.1 — 10 ao). 
In order to take the EFE into account, we must simply use 
the corresponding 7moat £>(??) function which now reads : 



jMONDiv) = 




>(B6) 



with g N {n) = GM{n)/rf. Comparing Eq.(B6) with (14) 
shows that, for given stellar parameters, ^mond (jf) for EFE 
is smaller than the one for the simple function so that the 
predicted velocity dispersion is smaller. Such a somewhat 
counterintuitive result can be qualitatively explained by not- 
ing that the EFE term g ext increases the Newtonian accel- 
eration and delays the transition from the Newtonian to the 
MOND regime. As a consequence, the boost in acceleration 
(and hence in ao) due to the MOND effect starts later and 
reduces ao with respect to the no EFE case. 

In order to study the impact of the EFE on the MFP co- 
efficients, we start from the model with constant anisotropy 
and no EFE and rewrite Eq.(27) as : 

log Reff = a log I EFE 

x a% FE ^j + b log (E e ) + c log n + d 

with ao and aQ FE the velocity dispersion without and with 
the EFE taken into account. Note that proceeding this way 
allows us to estimate (a, 6, c, d) using the values found for 
the no EFE simple case. In a first reasonably good approx- 
imation, it is : 



(TO 



dext log ( ] + b e 

V do / 



log 

7o J \ a 

with (a ex t,bext) — (0.09,0.21). If we further postulate that 
g e xt correlate with luminosity as 

log (gext/ao) = a ext logL + /3 ex t 

and repeat the same steps leading from Eq.(27) to Eq.(30), 
we find that the slope coefficients of the MFP now read : 



CtMFP — 



/3iV 



1 — 2(a*6 + a„c + aa ext OL ext ) 
(1 + a*)b + a v c + aa ext a ext 



(B7) 



1 — 2(a*b + a v c + aa ex ta ex t) 

which reduce to Eq.(31) for a ext — 0, i.e. no EFE effect. In 
order to get an estimate of a ext , we can assume that the 
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values of (a, b, c) are the same as those computed using the 
median /3 value obtained by fitting the constant anisotropy 
model with no EFE to the data. Setting (a*,av) to the 
median values in Table 1 for the i' filter and solving ckmfp = 
a b s , we finally get a ext ~ — f .5, i.e. brighter galaxies should 
be affected by a lower EFE. 

Investigating whether such a correlation is observation- 
ally motivated is a rather difficult task. In order to estimate 
g e xt, one can not resort to a fit to galaxies binned in lumi- 
nosity since the magnitude of the EFE depend on the envi- 
ronment in which a given galaxy is embedded. Although it is 
possible to study the environment of the ETGs in our sam- 
ple on a case -by -case basis, their large number makes this 
task quite beyond the scope of this paper. As a general re- 
mark, we however note that the EFE is typically invoked to 
improve the fit to the rotation curves of dwarf spiral galax- 
ies (see, e.g., Angus 2008), while the value of g ext is smaller 
for Milky Way like systems. This observation goes in the 
right direction. Moreover, one has to take care also of the 
symmetry of the problem. Brighter galaxies are typically 
more massive and hence populate the inner regions of clus- 
ters where they feel the combined action of many galaxies 
around. One can naively expect that the random orientation 
of the EFE from each companion leads to a sort of compen- 
sation thus lowering g ext for brighter systems. However, a 
further quantitative analysis is mandatory in order to lend 
support to this scenario. 
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